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'''deceased 

Multiplicity fluctuations of positively, negatively and all charged hadrons in the forward hemi- 
sphere were studied in central Pb-|-Pb collisions at 20^1, 30^, 40^1, 80^1 and 158^1 GeV. The mul- 
tiplicity distributions and their scaled variances lj are presented in dependence of collision energy 
as well as of rapidity and transverse momentum. The distributions have bell-like shape and their 
scaled variances are in the range from 0.8 to 1.2 without any significant structure in their energy 
dependence. No indication of the critical point in fluctuations are observed. The string-hadronic 
model UrQMD significantly overpredicts the mean, but approximately reproduces the scaled vari- 
ance of the multiplicity distributions. The predictions of the statistical hadron-resonance gas model 
obtained within the grand-canonical and canonical ensembles disagree with the measured scaled 
variances. The narrower than Poissonian multiplicity fluctuations measured in numerous cases may 
be explained by the impact of conservation laws on fluctuations in relativistic systems. 
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I. INTRODUCTION 

In matter of high energy densities (~ 1 GeV/fm'^) 
a phase transition is expected between hadrons and a 
state of quasi- free quarks and gluons, the quark gluon 
plasma (QGP) [H, [|. Measurements indicate that this 
critical energy density is exceeded at top SPS ^ 0] and 
RHIC d, H, H Hi energies during the early stage of heavy 
ion collisions. Moreover, the energy dependence of var- 
ious observables shows anomalies at low SPS energies 
which suggest the onset of deconfinement around 30^ 
GeV beam energy in central Pb-|-Pb collisions (ol. [lol. [lll| . 

It was predicted ^T5| that the onset of deconfinement 
can lead to a non-monotonic behaviour of multiplicity 
fluctuations. Lattice QCD calculations suggest further- 
more the existence of a critical point in the phase dia- 
gram of strongly interacting matter which separates the 
line of first order phase transition at high baryo-chemical 
potentials and low temperature from a crossover at low 
baryo-chemical potential and high temperature. An in- 
crease of multiplicity fluctuations near the critical point 
of strongly interacting matter is expected • 

In statistical models the widths of the multiplicity dis- 
tributions depend on the conservation laws which the 
system obeys. Even though for different statistical en- 
sembles the mean multiplicity is the same for sufficiently 
large volumes this is not necessarily so for higher mo- 
ments of the multiplicity distribution hence multiplicity 
fluctuations Fluctuations are largest in the grand- 
canonical ensemble, where all conservation laws are ful- 
filled only on average and not on an event-by-event ba- 
sis. The multiplicity fluctuations are much smaller in 
the canonical ensemble, where the electric and baryonic 
charges as well as strangeness are globally conserved. 
The smallest fluctuations are obtained within the micro- 
canonical ensemble, for which the charges as well as to- 
tal energy and momentum are conserved. It should be 
underlined that in non-relativistic gases the situation is 
very different, namely particle number is conserved in 
the micro-canonical and canonical ensembles and conse- 
quently the total multiplicity in these ensembles does not 
fluctuate. 

These theoretical considerations motivated vigorous 
theoretical [3i [Hi [Hi [l3i [3 ^^'^ experimental studies of 
multiplicity fluctuations in high energy nuclear collisions. 

Results on the centrality dependence of multiplic- 
ity fluctuations in Pb-|-Pb collisions obtained by the 
NA49 and WA98 ^ collaborations at top SPS en- 
ergy show an increase of multiplicity fluctuations with 
decreasing centrality of the collision in the forward hemi- 
sphere. A similar increase of multiplicity fluctuations is 
observed at midrapidity by the PHENIX [lH [13] collab- 
oration at RHIC energies. 

Transverse momentum fluctuations [23l | also show a 
non-monotonic dependence on system size. They in- 
crease from p+p to Si-|-Si and peripheral Pb-|-Pb col- 
lisions and decrease from peripheral to central Pb-|-Pb 
collisions. Possible relations to multiplicity fluctuations 



are discussed in [2J, |25| . Preliminary results of NA49 on 
the e nerg y dependence of transverse momentum fluctua- 
tions [2^ in central Pb-|-Pb collisions indicate a constant 
behaviour. 

This paper presents the dependence of multiplicity 
fluctuations on energy as well as on rapidity and trans- 
verse momentum for the most central Pb-|-Pb collisions 
at 20 A, 30 A, 40 A, 80 A and 158 A GeV as measured by 
the NA49 experiment at the CERN SPS. 

The paper is organized as follows. In chapter [TT] the 
notation and definitions are presented. In chapter IIIII 
the NA49 experiment and the experimental procedure of 
selecting events and tracks used for this analysis is de- 
scribed. In chapter llVI the experimental results on multi- 
plicity fluctuations are shown as a function of energy, ra- 
pidity and transverse momentum [2^ . These results are 
compared to the predictions of the hadron-resonance gas 
model and the string-hadronic model UrQMD [28[ 
in chapter |Vl Furthermore the measurements are also 
discussed with respect to the search for the onset of de- 
confinement and the critical point. The paper ends with 
a summary in chapter IVII 



II. MEASURE OF MULTIPLICITY 
FLUCTUATIONS 

Let P{n) denote the probability to observe a particle 
multiplicity n P{n) = 1) in a high energy nuclear 
collision. 

The scaled variance uj used in this paper as a measure 
of multiplicity fluctuations is commonly used in elemen- 
tary and heavy ion collisions, both for theoretical (see 
e.g. Ref s . fm [itI I29I [soj l and experimental (see e.g. 



Refs. [il, [lalJlHlfystudies. It is defined as 

_ Var{n) _ (n^) - (n)" 



(1) 



where Var(n) = X]n("'~ ^-^^^ ("-) = 12n '^^'("•) 

are variance and mean of the multiplicity distribution, 
respectively. 

In a superposition model uj is the same in A + A colli- 
sions as in nucleon-nucleon interactions at the same en- 
ergy per nucleon provided the number of particle pro- 
ducing sources does not fluctuate from event to event. 
String-hadronic models predict similar values of lo for 
p-l-p and Pb-I-Pb collisions |l3i (lO]- In a hadron-gas 
model [IB] the scaled variance converges quickly to a con- 
stant value with increasing volume of the system. In the 
special case of a hadron-gas model in the grand-canonical 
formulation jlQ] . neglecting quantum effects and reso- 
nance decays, the multiplicity distribution is a Poisson 
one, namely 



P{n) 



-in) 



(2) 



The variance of a Poisson distribution is equal to its 
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mean, and thus the scaled variance is w = 1, indepen- 
dent of mean multiplicity. 

If there are no particle correlations in momentum space 
and the single particle distribution is independent of par- 
ticle multiplicity the scaled variance of an arbitrary mul- 
tiplicity distribution observed in a limited acceptance 
is related to the scaled variance in the full phase-space 
("47r") as (see appendix [XT] and Refs. pj,, Jfi] for deriva- 
tion): 



(tJ47r ~l)p+l, 



(3) 



where p denotes the fraction of particles measured in the 
corresponding acceptance. Note that the dependence de- 
scribed by Eq. [3] is violated if effects like resonance de- 
cays, quantum statistics and energy- momentum conser- 
vation introduce correlations in momentum space [32]. 

In the following the scaled variances of the multiplic- 
ity distributions of positively, negatively and all charged 
hadrons are denoted as uj{h~^), u)[h~) and respec- 
tively. 



III. THE NA49 EXPERIMENT 

The NA49 detector [sl] (see Fig. [T]) is a large accep- 
tance fixed target hadron spectrometer. Its main devices 
are four large volume time projection chambers (TPCs). 
Two of them, called vertex-TPCs (VTPC-1 and 2), are 
located in two superconducting dipole magnets (VTX-1 
and 2) with a total bending power up to 7.8 Tm. The 
magnetic field used at 158A GeV (B(VTX-l)w 1.5 T) 
and B(VTX-2)w 1.1 T) was scaled down in proportion 
to the beam energy for lower energies. The other two 
TPCs (MTPC-L and MTPC-R), called main-TPCs, are 
installed behind the magnets on the left and the right side 
of the beam line allowing precise particle tracking. The 
measurement of the energy loss dE/dx in the detector 
gas provides particle identification in a large momentum 
range. It is complemented by time of flight (TOF) detec- 
tors measuring particles at mid-rapidity. In this analysis 
dE/dx information is used only to reject electrons. 

The target is located 80 cm upstream of the first vertex 
TPC. The target thickness is 0.2 mm (0.224 g/cm^) for 
2QA - 80A GeV and 0.3 mm (0.336 g/cm^) for 158A 
GeV. Using 7.15 barn as the inelastic cross-section for 
Pb-I-Pb collisions this yields an interaction probability of 
0.46% and 0.7%, respectively. The interaction length of 
the strong interaction for Pb ions in a Pb target is 4.26 
cm. 

Three beam-position-detectors (BPDs) allow a precise 
determination of the point where the beam hits the target 
foil. The centrality of a collision is determined by measur- 
ing the energy of projectile spectators in the downstream 
veto calorimeter (VCAL, see section UlI B[) . The accep- 
tance of the veto calorimeter is adjusted at each energy 
by a proper setup of the collimator (COLL). 



energy (GoV) 


number of events 


2QA 


6602 


30^ 


8219 


40^ 


21995 


80^1 


2307 


158yl 


5493 



TABLE I: Statistics for the 1% most central collisions used 
for this analysis at different beam energies. 



A. Data Sets and Event Selection 

In this publication the results for central Pb-f Pb col- 
lisions at 20v4, mA, AOA, 80A and 158A GeV are pre- 
sented. The numbers of events used from these data sets 
are given in Table |T1 

In order to get a "clean" sample of events excluding 
for instance collisions outside the target or event pileup, 
the following event selection criteria are applied to data: 

• The fit of the interaction point, based on the recon- 
structed tracks, was successful. 

• The position of the fitted interaction point is close 
to the position obtained from the beam position 
detectors. 

• At least 10% of all tracks are used for the recon- 
struction of the interaction point. The reconstruc- 
tion of the interaction point was optimized for pre- 
cision by selecting long and well measured tracks 
in an iterative procedure. 

The event cuts have a small influence on uj, the results 
differ by less than 1% when only the cut requirement of 
a successful fit of the main vertex is used. 

Beam lead ions which do not interact strongly in the 
target produce delta electrons both in the target foil and 
the detector gas. These electrons may curl up in the 
TPCs, increase their occupancy and might therefore re- 
duce the reconstruction efficiency. In order to avoid this 
effect only those events are selected for the analysis in 
which there are no beam ions passing through the detec- 
tor within the read-out time of the event. 



B. Centrality Selection 

Fluctuations in the number of participants lead to an 
increase of multiplicity fluctuations. In a superposition 
model the total multiplicity n is the sum of the number 
of particles produced by k particle production sources: 



(4) 



where the summation index i runs over the sources. Un- 
der the assumption of statistically identical sources the 
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FIG. 1: (Color online) Setup of the NA49 experiment for Pb+Pb collisions, see text for more details. 
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FIG. 2: A sketch of the horizontal deflection for charged parti- 
cles at the front face of the iron collimator for the 158^ GeV 
magnetic field setting. The broadened distribution of each 
species is due to the Fermi motion of nucleons or fragments; 
additionally, the oval shapes are due to the deflection of 
charged particles in the magnetic field. The sizes of the dis- 
tributions correspond to one standard deviation. The open 
circles in the fragment acceptance represent particles of Z/A 
other than one half [3^ . 

scaled variance lo of the multiplicity distribution has two 
contributions. The first is due to the fluctuations of the 
number of particles emitted by a single source uJso, the 
second is due to the fluctuations in the number of sources 
ujk (see appendix lA 21 for derivation): 

uj^oj-"' + {n"") ■ ujk, (5) 

where (n'*°) is the mean multipHcity of hadrons from a 
single source. The fluctuations in the number of sources 
ujk can be attributed to fluctuations in the number of pro- 
jectile and target participants. In order to minimize the 
fluctuations of the number of participants the centrality 
variation in the ensemble of events should be as smaU as 
possible, for which very central collisions are best suited. 

In order to fix the number of projectile participants 
the NA49 experiment uses the energy in the projectile 
spectator domain as a measure of centrality, called pro- 



energy 


collim 
X (cm) 


ator 
y (cm) 


ring calorimeter 
X (cm) 


20 








10 


30 








10 


40 


-13 


+47 


±12 


17 


80 


-13 


+47 


±12 


17 


158 


-5 


+38 


±5 


17 



TABLE II: Settings of the collimator and the ring calorimeter 
defining the acceptance of the veto calorimeter for different 
energies with respect to the position of neutrons with zero 
transverse momentum. See the text for more details. 



jectile centrality below. The downstream veto calorime- 
ter [3^ of NA49, originally designed for NA5, measures 
the energy carried by the particles in the projectile spec- 
tator phase space region |3^ . A collimator in front of the 
calorimeter is located 25 m downstream from the target 
and is adjusted for each energy in such a way that all 
projectile spectator protons, neutrons and fragments can 
reach the veto calorimeter. For 158A GeV the hole in 
the collimator extends ±5 cm in vertical direction and 
—5 cm and +38 cm in horizontal direction taking into 
account the deflection of charged spectators by the mag- 
netic field (Fig. [21 Table [n|. Due to a larger spread of 
spectators, the hole of the collimator is larger for 40^4 and 
80 A GeV. For 20 A and 30 A the collimator is removed and 
the ring calorimeter (RCAL in Fig. [1]) positioned 18 m 
downstream from the target serves as a collimator. 

The settings of the hole in the collimator and the po- 
sition of the ring calorimeter for the different energies is 
shown in Table |lll The zero point is the point where 
neutrons with no transverse momentum would pass the 
collimator. The collimator is not symmetric around the 
zero point because the nuclear fragments and spectator 
protons carry positive charge and are deflected by the 
magnetic fleld in positive x direction. The last column in 
the table is the position of the center of the ring calorime- 
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ter. Its hole has a radius of 28 cm. 

The acceptance of the veto calorimeter for neutral and 
positive particles for 158^ GeV is shown in Fig. [51 A c- 
ceptance tables in p, pT and (p can be obtained at [36] . 

Due to the geometry of the collimator and the mag- 
netic field, a small number of positive and neutral non- 
spectator particles can hit the veto calorimeter. For 
positively charged particles, the acceptance of the TPCs 
and the veto calorimeter overlap partly. The maximum 
amount of a possible auto-correlation is estimated by a 
comparison of Ci;(/i+) for UrQMD events selected by their 
veto energy to UrQMD events with a zero impact pa- 
rameter in the forward region (Fig. I12p and found to be 
smaller than 3%. 

The acceptance of the veto calorimeter for negatively 
charged particles is very small because they are bent by 
the magnetic field into the direction opposite to the one 
of the positively charged particles, and the collimator is 
adjusted to detect positively charged and neutral projec- 
tile spectators. 

The projectile centrality Cproj of an event with a veto 
energy Eveto is defined as the percentage of all inelas- 
tic events which are as central or more central than the 
given event according to the energy deposited in the veto 
calorimeter by the projectile spectator nucleons. Smaller 
Cproj correspond to more central events. Using the frac- 
tion of inelastic cross section Ctrig = accepted by 
the trigger {atrig is derived from the target thickness and 
the interaction rate, amei is assumed to be 7.15 barn) and 
the veto energy distribution C'proj is given by: 



C 



Proj 



trig 



Jo dN /dEveto,trigdEveto 



(6) 



where dN/dEyeto,trig is the veto calorimeter energy dis- 
tribution for a given trigger. 

The finite resolution of the veto calorimeter causes 
additional fluctuations in the number of participants. 
Based on the analysis of the NA49 Pb-|-Pb data the res- 
olution of the veto calorimeter was estimated in [l^ to 
be: 



a{E- 



Veto) 



2.85 



16 



Eve 



VE- 



Veto 



Eve 



(7) 



where Eveto is in units of GeV. In order to check this 
parametrization, the distribution of the spectators was 
simulated by the SHIELD model [s;^. The SHIELD 
model delivers both spectator nucleons and nuclear frag- 
ments, in contrast to most string hadronic models, which 
only produce spectator nucleons. A simulation per- 
formed at 20 A and 158 A GeV including the geometry 
of the NA49 detector and the non-uniformity of the veto 
calorimeter confirms the parametrization given by Eq. [7| 
as an upper hmit (see Fig. [4]). 

The veto calorimeter response can in principle change 
with time (aging effects, etc.). Therefore a time depen- 
dent calibration of the veto energy was applied. The 
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FIG. 3: Acceptance of the veto calorimeter for neutral (top) 
and positively charged (bottom) main vertex particles at 
158A GeV as a function of total momentum p and transverse 
momentum pr- 



contribution of this correction to uj turned out to be very 
smaU (< 1%, see Table |lVl. 

When fixing the projectile centrality Cproj (Eq- [6]), 
thereby fixing the number of projectile participants 



N 



Proj 



, the number of target participants Np""^^ can still 
fluctuate. Thus the total number of participants is not 
rigorously constant and could contribute to fluctuations. 
The fluctuations of the number of target participants ob- 
tained by UrQMD and HSD simulations [sHi, expressed 

as their scaled variance 



Targ 



rTarg 



are shown in Fig. [5j For non-central collisions the num- 
ber of target participants strongly fluctuates, even for a 
flxed number of projectile participants. This is con- 
sistent with the increase of lj with decreasing centrality 
observed in the forward hemisphe re [l9l [39| . However, 
alternative explanations also exist [2^, Eo] ■ 

For further analysis the 1% most central collisions (ac- 
cording to their veto energy) are selected in order to min- 
imize the fiucutations in the number of participants. For 
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Eve.o(GeV) 

FIG. 4: Resolution of the veto calorimeter estimated 
by a SHIELD simulation (histogram) compared to the 
parametrization Eq.[7](solid line) for 20^4 (top) and 158^4 GeV 
(bottom). 



these very central collisions, the fluctuation in the num- 
ber of target participants is expected to be smallest and 
its scaled variance luJ,""^^ is expected to be about 0.1 (see 
Fig. [5]) for an estimated number of target participants of 
iVj'"'^ « 192. 

In order to estimate the effect on uj of target partici- 
pant fluctuations and non-spectator particles in the veto 
calorimeter, the energy dependence of the scaled vari- 
ance of the multiplicity distribution is calculated in the 
UrQMD 1.3 model both for collisions with zero impact 
parameter and for collisions selected according to their 
veto energy. The resulting difference of lo in the forward 
acceptance (see section UlI CI) is smaller than 2% for neg- 
atively, smaller than 3% for positively and smaller than 
4% for all charged hadrons. In the midrapidity region the 
influence of the fluctuations of target participants on to is 
expected to be much larger. Indeed, the differences of uj 
increase to up to 6% for negative, up to 9% for positive 
and up to 13% for all charged hadrons. 




FIG. 5: (Color online) Scaled variance of the number of target 
participants for a fixed number of projectile participants in 
the UrQMD and HSD models. The plot is taken from 3]. 

In order to check the influence of the centrality selec- 
tion, UJ was also determined for the 0.5% most central 
collisions. The change compared to the values obtained 
for the 1% most central collisions is smaller than 3% for 
positive, 2% for negative and 5% for all charged hadrons. 



C. Track Selection 

Since detector effects like track reconstruction effi- 
ciency might have a significant infiuence on multiplicity 
fluctuations, it is important to select a sample of well 
defined tracks for the analysis. The following track selec- 
tion criteria are used for this analysis and are explained 
in this section: 

• Number of potential points (the number of points 
a track can have according to its geometry) in the 
TPCs: > 30. 

• The ratio of the number of reconstructed points to 
the number of potential points: > 0.5. 

• Sum of the number of reconstructed points in 
VTPC-1 and 2: > 5. 

• Sum of the number of reconstructed points in 
VTPC-2 and MTPCs: > 5. 

• The track is extrapolated to the plane of the target 
foil. This point must be closer than 4 cm in x- and 
2 cm in y- direction to the interaction point of the 
collision. 
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• In order to exclude electrons from the analysis, a 
cut on the energy loss [dE/dx) in the detector gas 
was applied. All tracks with an energy loss more 
than 0.2 minimum ionising units higher than the 
pion dE/dx (in the region of the relativistic rise of 
the Bethe-Bloch formula) are rejected. 

The reconstruction efficiency is calculated using the em- 
bedding method. Events containing a few tracks were 
generated and processed by the simulation software. The 
resulting raw data were embedded into real events. The 
combined raw data were reconstructed and the input 
tracks were matched with the reconstructed ones. Em- 
bedding simulations show a significant decrease of recon- 
struction efficiency with increasing event multiplicity in 
the midrapidity region at 158A GeV using the track selec- 
tion criteria described above. Therefore for this energy 
an additional cut was used, namely that tracks should 
have at least 5 reconstructed points both in VTPC-2 and 
in the MTPCs. For these tracks no significant depen- 
dence of reconstruction efficiency on track multiplicity is 
observed. 

Reconstruction inefficiencies mostly occur for tracks 
with a very low number of points in the TPCs or for 
tracks which only have points in the VTPC-1 or in the 
main TPC. These tracks are not used for this analysis. 

In the following the longitudinal motion of particles 
is characterized by the rapidity in the center of mass 
system assuming pion mass of the particle. This measure 
is called pion rapidity and is denoted as 2/(7r). 

The distributions of the registered tracks after apply- 
ing the track selection criteria are shown in Fig. [S] as a 
function of pion rapidity 2/(7r) and transverse momentum 
Pt- Acceptance tables in yin), pt and (j) can be obtained 
from [sHi . Only tracks in the rapidity interval starting at 
midrapidity and ending at beam rapidity are used. 

In order to study the multiplicity fluctuations differ- 
entially, the pion rapidity interval < j/(7r) < yteam is 
divided into two parts, the "midrapidity" (0 < y{TT) < 1) 
and the "forward rapidity" (1 < y{Tr) < ybeam) region 
(see Fig. [71). The fractions of total charged particle mul- 
tiplicity falling into the different rapidity intervals are 
given in Table IIIII and Fig. [51 The values are calculated 
using the VENUS event generator [4lj as input for a 
GEANT simulation. The tracks produced by GEANT 
are converted into detector signals and reconstructed by 
the NA49 reconstruction chain. For the determination 
of the acceptance the negatively charged main vertex pi- 
ous, kaons and anti-protons are used. In both regions 
a similar number of particles is detected by NA49. In 
the forward acceptance the particles are mostly passing 
through both the vertex- and the main- TPCs and are 
therefore efficiently reconstructed for all collision ener- 
gies. According to the UrQMD model the fluctuations 
in the number of target participants contribute mostly 
to the particle number fluctuations in the target hemi- 
sphere and the midrapidity region. Their influence on 
uj in the forward region (j/(7r) > 1) can be estimated by 
the difference in scaled variance between 6 = and veto 
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FIG. 6: Distribution of detected negatively charged particles 
which fulfill the track selection criteria as a function of y(7r) 
and PT for 2QA (top), 30^, 40A, 80.4 and 158A GeV (bot- 
tom). 
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FIG. 7: (Color online) Dashed line: Double-Gauss 
parametrization of the rapidity distribution of negatively 
charged pions and kaons in Pb+Pb collisions at 20^4 (top) [10] 
and 158^4 GeV (bottom) '9]. The solid line is the measured 
y(7r) distribution with the track selection criteria described 
in section IIII CI The vertical lines indicate the limits of the 
rapidity intervals y{Tv) = 0, y{-!v) = 1 and j/(7r) = ybeam used 
for this analysis. 
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FIG. 8: (Color online) Fraction p of total negatively charged 
main vertex pion, kaon and anti-proton multiplicity which is 
accepted and reconstructed as a function of collision energy. 
Circles: < y{Tr) < ybeam, boxes: < y{n) < 1, triangles: 

1 < y{Tv) < ytearn- 



energy 


< y{n) 

ybeam 


< y{n) 
< 1 


1 < y{n) 

ybeam 


a{y){TT ) 


20 


15.3% 


7.2% 


8.1% 


1.01 


30 


19.1% 


8.4% 


10.7% 


1.08 


40 


21.7% 


9.2% 


12.6% 


1.1 


80 


28.2% 


11.2% 


17% 


1.23 


158 


28.8% 


9.6% 


19.2% 


1.38 



TABLE III: Fraction (in percent) of negatively charged main 
vertex pions, kaons and anti-protons in different rapidity in- 
tervals for different collision energies which are accepted and 
reconstructed. In addition, the width of the rapidity distri- 
bution of negatively charged pions is given [gl. [lo|. 



selected collisions (see section UlI Bp and is about 1 — 2%. 

Note that the acceptance used for this analysis is larger 
than the one used for the preliminary data shown in [42l . 




D. Systematic Errors 

The influence of the selection criteria described above 
on the scaled variance lu of the multiplicity distribution 
has been studied and the results are presented in Table HVl 
and Figs. [MTT] The event selection criteria described in 
section lTlI Al change uj by up to 2% compared to the value 
obtained when not applying these cuts. The finite resolu- 
tion of the veto calorimeter causes additional fluctuations 
in the number of projectile participants and therefore in- 
creases the measured lo. In a superposition model the 
effect of the veto calorimeter resolution is estimated to 



be [3: 

where Efy^am is the total energy per projectile nuclcon. 
The parametrization Eq. ([7]), which serves as an upper 
limit of the resolution of the calorimeter, was used to de- 
termine the potential influence of the resolution on lo. 
For the very central collisions selected for this analy- 
sis the measured uo is found to increase due to the fi- 
nite calorimeter resolution by less than 1.5%. Therefore 
a correction for this effect is not applied. In order to 
take possible aging effects of the calorimeter (see sec- 
tion |TTTB]) into account, a time dependent calibration is 
applied to the measured veto energy. However, the ef- 
fect of this correction is very small, uo changes by less 
than 1%. Track selection criteria are applied to remove 
electrons and tracks not originating from the main in- 
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TABLE IV: Maximum change Alu of the scaled variance uj of 
the multiplicity distribution for positively, negatively and all 
charged hadrons when applying a correction or neglecting a 
cut. The systematic errors are calculated by adding the error 
contributions in quadrature. The last row shows the change 
of CO resulting from a change in the centrality selection from 
1% to 0.5%. 



ter action point. The value of u) is changed by less than 
1.5% for positively and negatively and less than 3% for 
all charged hadrons when removing these cuts. 

Embedding simulations demonstrated that the recon- 
struction efficiency shows no significant decrease with in- 
creasing particle multiplicity. Therefore no systematic er- 
ror due to reconstruction efficiency was attributed. The 
overall reconstruction efficiency is about 95% and is in- 
cluded in the calculation of the acceptances (Fig. [51 Ta- 
ble |TTl|. 

The total systematic error is calculated by adding the 
contributions of the different error sources in quadrature. 
It is 2.4%, 1.8% and 3.8% for positively, negatively and 
all charged hadrons, respectively. 

In order to estimate the effect of centrality selection, 
also the 0.5% most central collisions are studied. The 
result for ui for this stricter selection is up to 5% different 
from that obtained for the 1% most central collisions. 
As the centrality selection is a well defined procedure 
and can be repeated in model calculations, the difference 
of uj for the 0.5% and 1% most central collisions is not 
considered as part of the systematic error. 



IV. RESULTS ON MULTIPLICITY 
FLUCTUATIONS 
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FIG. 9: (Color online) Systematic errors and difference be- 
tween the 0.5% and the 1% most central collisions of the 
scaled variance uj of the multiplicity distribution for positively 
charged hadrons at midrapidity (0 < J/(7r) < 1, top) and for- 
ward acceptance (1 < y(7i") < ybeam, bottom) as a function 
of collision energy. Lu{std.) corresponds to the value obtained 
when using the standard event and track selection criteria and 
no correction for the veto calorimeter resolution. 



In this chapter results on multiplicity fluctuations for 
negatively, positively and all charged hadrons are pre- 
sented for Pb-I-Pb colhsions at 20A, AOA, 80A and 
158A GeV. In order to minimize the fluctuations in the 
number of participants, the 1% most central collisions ac- 
cording to the energy of projectile spectators measured in 
the veto calorimeter are selected (see section UlI B[) . The 
rapidity interval < y(7r) < ybeam used for this analysis 
is divided into two subintervals, < y(7r) < 1 ("midra- 
pidity") and 1 < y(7r) < ybeam ("forward rapidity", see 
section ImC]) . 

In the following the errors indicated by vertical lines 
with attached horizontal bars correspond to the statisti- 
cal errors only, the thick horizontal bars are the statistical 



and systematic errors added in quadrature. 



A. Multiplicity Distributions 

The multiplicity distributions for the different ener- 
gies, charges and rapidity intervals as well as the ratios 
of the measured multiplicity distributions to a Poisson 
distribution with the same mean multiplicity are shown 
in Figs. [26ll34l For the ratio to the Poisson distribu- 
tions only points with statistical errors smaller than 20% 
are shown. All multiplicity distributions have a bell- 
like shape, and no significant tails or events with a very 
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FIG. 10: (Color online) Systematic errors and difference be- 
tween the 0.5% and the 1% most central collisions of the 
scaled variance uj of the multiplicity distribution for nega- 
tively charged hadrons at midrapidity (0 < 3/('7r) < 1, top) 
and forward acceptance (1 < y{iT) < ybeam, bottom) as a 
function of collision energy. uj{std.) corresponds to the value 
obtained when using the standard event and track selection 
criteria and no correction for the veto calorimeter resolution. 
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FIG. 11: (Color online) Systematic errors and difference be- 
tween the 0.5% and the 1% most central collisions of the scaled 
variance to of the multiplicity distribution for all charged 
hadrons at midrapidity (0 < y{Tr) < 1, top) and forward 
acceptance (1 < j/(7r) < yteam, bottom) as a function of colli- 
sion energy. uj(std.) corresponds to the value obtained when 
using the standard event and track selection criteria and no 
correction for the veto calorimeter resolution. d 



high or very low multiplicity are observed. The ratios of 
measured multiplicity distributions to the corresponding 
Poisson distributions are symmetric around their mean 
value. 

The measured multiplicity distributions are narrower 
than the Poisson ones in the forward acceptance for pos- 
itively and negatively charged hadrons at all energies. In 
the midrapidity acceptance the measured distributions 
are wider or similar to the Poisson ones. The distribu- 
tions for all charged hadrons are broader than the ones 
for positively and negatively charged particles separately. 



B. Energy Dependence of u 

The energy dependence of the scaled variance uj of the 
multiplicity distributions for negatively, positively and all 
charged particles for three rapidity intervals is shown in 
Figs. [T2lfT4l the numerical values are given in Table fVl 

For positively and negatively charged hadrons the val- 
ues of oj are similar and smaller than 1 in the very forward 
region (1 < y(7r) < ybeam) at all energies. At midrapid- 
ity they are larger than 1. For all charged particles lo is 
larger than for each charge separately. 

No significant structure or non-monotonic behaviour is 
observed in the energy dependence of uj. 
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FIG. 12: (Color online) Scaled variance uj of the multiplicity 
distribution of positively charged hadrons produced in cen- 
tral Pb+Pb collisions as a function of collision energy. Top: 
full experimental acceptance, middle: midrapidity, bottom: 
forward rapidity. 



FIG. 13: (Color online) Scaled variance uj of the multiplicity 
distribution of negatively charged hadrons produced in cen- 
tral Pb-(-Pb collisions as a function of collision energy. Top: 
full experimental acceptance, middle: midrapidity, bottom: 
forward rapidity. 



Signatures of the critical point are exp ected to oc- 
cur mostly at low transverse momenta [13|. The energy 
dependence of multiplicity fluctuations for low trans- 
verse momentum particles is shown in Fig. 1151 No non- 
monotonic behaviour is observed. 



C. Rapidity Dependence of uj 

The rapidity dependence of the scaled variance uj of 
the multiphcity distributions for 20 A, SO A, 40 A, 80 A and 
158A GeV central Pb-|-Pb collisions is shown in Figs.fT^ 
[T51 In order to remove the "trivial" dependence of u on 
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FIG. 14: (Color online) Scaled variance uj of the multiplic- 
ity distribution of all charged hadrons produced in central 
Pb+Pb collisions as a function of collision energy. Top: full 
experimental acceptance, middle: midrapidity, bottom: for- 
ward rapidity. 



the fraction of accepted tracks (see Eq. [3]) the rapidity 
bins j/c — Ay < y < jjc + Ay are constructed in such a 
way that the mean multiphcity in each bin is the same. 

If there were no correlations in momentum space and 
the single particle spectra are independent of particle 
multiplicity, the resulting values of lu shown in Figs. 1161 



FIG. 15: (Color online) Scaled variance ui of the multiplicity 
distribution of negatively charged hadrons with low trans- 
verse momentum at forward rapidities produced in central 
Pb-|-Pb collisions as a function of collision energy. Top: 
Pt < 0.3 GeV/c, bottom: pr < 0.5 GeV/c. 



1181 would be independent of rapidity. This is not the 
case, the experimental data show an increase of w to- 
wards midrapidity for all charges and energies. 



D. Transverse Momentum Dependence of lu 

The transverse momentum dependence of u at top SPS 
energy is shown in Fig. 1191 The transverse momentum 
range of — 1.5 GeV/c is divided into five bins in such a 
way that the mean multiplicity in each bin is the same. 
The horizontal position of the points in Fig. [19] corre- 
spond to the center of gravity of the transverse momen- 
tum distribution in the transverse momentum range of 
the corresponding bin. Only a small rapidity interval in 
the forward acceptance (1.25 < y(7r) < 1.75) is used for 
this study. A larger rapidity interval might cause a bias 
because the acceptance in rapidity is different for differ- 
ent transverse momenta. 

An increase of lu with decreasing transverse momen- 
tum, which is more pronounced for Lu{h~) than for aj(/i+). 



13 




1.2 



■ NA49 

^ • UrQMD 


1 T - % , 


- — i — 





0.5 1 1.5 2 




FIG. 16: (Color online) Rapidity dependence of the scaled 
variance u) of the multiplicity distribution of positively 
charged hadrons in central Pb+Pb collisions at 20^4 (top), 
30A, 40A, 80A and 158A GeV (bottom) compared to UrQMD 
predictions with a centrality selection similar to the one for 
the experimental data. The rapidity bins are constructed in 
such a way that the mean multiplicity in each bin is the same. 
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FIG. 17: (Color online) Rapidity dependence of the scaled 
variance lo of the multiplicity distribution of negatively 
charged hadrons in central Pb+Pb collisions at 20 A (top), 
30A, 40A, 80A and 158A GeV (bottom) compared to UrQMD 
predictions with a centrality selection similar to the one for 
the experimental data. The rapidity bins are constructed in 
such a way that the mean multiplicity in each bin is the same. 
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FIG. 19: (Color online) Transverse momentum dependence of 
the scaled variance of the multiplicity distribution of positive 
(top), negative and all charged (bottom) hadrons in the ra- 
pidity interval 1.25 < j/(7r) < 1.75 in central Pb+Pb collisions 
at 158A GeV . 
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FIG. 18: (Color online) Rapidity dependence of the scaled 
variance to of the multiplicity distribution of all charged 
hadrons in central Pb+Pb collisions at 20^1 (top), 30 A, 40 A, 
80A and 158A GeV (bottom) compared to UrQMD predic- 
tions with a centrality selection similar to the one for the 
experimental data. The rapidity bins are constructed in such 
a way that the mean multiplicity in each bin is the same. 



is found. Only the top SPS energy is shown because at 
lower energies the azimuthal acceptance of the NA49 de- 
tector is much smaller and therefore lu would approach 
one due to the small multiplicity. 
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FIG. 20: (Color online) Predictions of a hadron-resonance gas 
model for the scaled variance oj of the multiplicity distribution 
in full phase space for negatively charged hadrons. The pa- 
rameters of the ensemble, (T, /is) are the values of the chemi- 
cal freeze-out obtained by a hadron-gas model fit to produced 
particle ratios at different energies. Results are shown for the 
grand-canonical (GCE), canonical (CE) and micro-canonical 
ensemble (MCE). The plot is taken from Ref. 



V. MODEL COMPARISON 
A. Hadron-resonance gas model 

In a hadron-resonance gas model an equilibrium state 
of hadrons and hadronic resonances is assumed. Three 
different statistical ensembles are considered, namely the 
grand-canonical, the canonical and the micro-canonical 
ensemble, which differ by the conservation laws which 
are taken into account. In the grand-canonical ensemble 
conservation laws are not obeyed on an event-by-event 
basis, whereas in the canonical ensemble the total baryon 
number, strangeness and electrical charge have to be con- 
served in each event. In the micro-canonical ensemble the 
total energy and momentum are conserved in addition. 

In [l^ the fluctuations of particle multiplicity in full 
phase-space were calculated for these three different sta- 
tistical ensembles in the infinite volume limit. The energy 
dependence of multiplicity fluctuations is introduced via 
the chemical freeze-out parameters T (temperature) and 
/is (baryo-chemical potential), which have been deter- 
mined by hadron-resonance gas model fits at all energies 
to the mean particle multiplicities. Quantum statistics 
and resonance decays are included in the model calcula- 
tions. The scaled variance oj of the multiplicity distribu- 
tion of negatively charged hadrons is shown in Fig. as 
a function of collision energy. 

The results for oj in the micro-canonical, canonical and 
grand canonical ensemble are very different at high col- 
lision energies. The well known equivalence of statistical 
ensembles in the large volume limit only holds for mean 



values, not for multiplicity fluctuations. 

The value of lo is the largest in the grand-canonical 
ensemble. In the micro-canonical ensemble it is the 
smallest, the canonical ensemble lies in between. In 
the canonical and micro-canonical ensemble for posi- 
tively and negatively charged particles separately nar- 
rower than Poisson {lo < 1) multiplicity fluctuations are 
expected. The difference between the grand-canonical, 
canonical and micro-canonical ensemble show the impor- 
tance of a proper treatment of conservation laws for mod- 
elling multiplicity fluctuations. 

In order to compare the hadron resonance gas model 
predictions with experimental data, u) calculated in full 
phase space is extrapolated to the experimental accep- 
tance using Eq. [31 Although quantum effects and reso- 
nance decays introduce correlations in momentum space, 
Eq.[3]is the only presently known way to compare the pre- 
dictions of the grand-canonical and canonical ensemble to 
the experimental data. For the micro-canonical ensem- 
ble the energy and momentum conservation introduces 
stronger correlations in momentum space {s^. There- 
fore Eq. [3] cannot serve as a reasonable approximation. 
Resonance decays introduce only a weak correlation in 
momentum space for positively and negatively charged 
hadrons, because only a small number of resonances de- 
cay into two particles with the same charge. In con- 
trast a large number of resonances decay into two oppo- 
sitely charged hadrons, therefore Eq.[2]is not valid for all 
charged hadrons. 

At forward rapidity (1 < y{'K) < ybeam,', Figs. \T[\ and 
mi bottom), the fluctuations are overpredicted by both 
the canonical and the grand canonical models. However, 
the canonical model is closer to data. A micro-canonical 
ensemble predicts smaller fluctuations than the canonical 
model, but a quantitative comparison with data is not 
possible yet, because correlations in momentum space do 
not allow to extrapolate to the experimental acceptance 
using Eq. [31 

At midrapidity uj of the data (squares in Figs. ?IT\ and 
1221 top) is higher than in the forward region. In contrast 
to the experimental data the fluctuations in the number 
of target participants are not included in the hadron-gas 
model. From comparison of UrQMD simulations for b = 
collisions and collisions selected according to their veto 
energy it can be estimated that the target participant 
fluctuations increase lo by up to 9% in the midrapidity 
region. 

The shape of the measured multiplicity distribution is 
compared to the hadron-resonance gas model prediction 
for negatively charged hadrons at 158 A GeV in the for- 
ward acceptance in Fig. [23] (top). For this comparison 
the multiplicity distributions for the data and the model 
predictions are divided by Poisson distributions with the 
same mean multiplicities. The hadron-resonance gas 
model predicts a Gaussian-like shaped multiplicity distri- 
bution in full phase space [ij . Since this model gives no 
prediction about the mean multiplicity, it is taken from 
data. In order to calculate the multiplicity distribution 
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FIG. 21: (Color online) Scaled variance uj of the multiplicity 
distribution of positively charged hadrons produced in central 
Pb+Pb collisions as a function of collision energy in midra- 
pidity (top) and forward (bottom) acceptance compared to 
predictions of a grand canonical and canonical ensemble [l^ . 



FIG. 22: (Color online) Scaled variance uj of the multiplicity 
distribution of negatively charged hadrons produced in central 
Pb+Pb collisions as a function of collision energy in midra- 
pidity (top) and forward (bottom) acceptance compared to 
predictions of a grand canonical and canonical ensemble [igI ]. 



in the limited experimental acceptance the distribution 
in the full phase space is folded with a Binomial dis- 
tribution accepting the same fraction p of tracks as the 
experimental acceptance: 

[N — njW. 

where N is the multiplicity in the full phase space and 
n the multiplicity in the experimental acceptance. The 
multiplicity distribution in the experimental acceptance 
is given by 

) =^P4.(VV)BArN. (10) 

N 

Note that this procedure assumes that there are no cor- 
relations in momentum space. 

The ratio for the grand-canonical ensemble has a con- 
cave shape, i.e. the multiplicity distribution is wider than 
a Poisson distribution. For the canonical ensemble the 
shape is convex, showing that the distribution is nar- 



rower. The shape for the experimental data is more con- 
vex, demonstrating that the measured multiplicity dis- 
tribution is even narrower than the canonical one. 

In the canonical and grand canonical ensembles of the 
hadron-resonance gas model no mechanisms are present 
which would introduce a strong dependence of multi- 
plicity fluctuations on rapidity or transverse momentum, 
which is observed in the data and in UrQMD (Figs. [TB]- 
[T^ . In a three-pion gas statistical model using the 
micro-canonical ensemble an increase of fluctuations near 
midrapidity and for low pT was observed [s^l as an effect 
of energy and momentum conservation. 

B. String-hadronic models 

In this section the experimental data on mul- 
tiplicity fluctuations are compared to the outcome 
of string-hadronic model calculations, namely of the 
Ultra-relativistic Quantum Molecular Dynamics model 
(UrQMD vl.3) [H,!!!! and the Hadron-String Dynamics 
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FIG. 23: (Color online) Ratio of multiplicity distribution of 
experimental data (top) and UrQMD simulation (bottom) to 
a Poisson distribution with the same mean value for nega- 
tively charged hadrons in Pb+Pb collisions at 158A GeV in 
the forward acceptance. Only points with statistical errors 
smaller than 20% are shown. Hadron gas model predictions 
in the grand-canonical and canonical ensemble with the same 
mean multiplicity and fraction of accepted tracks are shown 
by lines. 



model (HSD) [465. 

The UrQMD microscopic transport approach is based 
on the propagation of constituent quarks and di-quarks 
accompanied by mesonic and baryonic degrees of free- 
dom. It simulates multiple interactions of in-going and 
newly produced particles, the excitation and fragmen- 
tation of colour strings and the formation and decay of 
hadronic resonances. Towards higher energies, the treat- 
ment of sub-hadronic degrees of freedom is of major im- 
portance. A phase transition to a quark-gluon state is 
not incorporated explicitly into the model dynamics. 

The scaled variance uj of the multiplicity distribution of 
negatively charged hadrons for all inelastic p+p and p-|-n 
interactions as well as central (6 = 0) Pb-|-Pb collisions 
predicted by the UrQMD model [l3| is shown in Fig. [21] 
in dependence of the collision energy. 

The scaled variance of multiplicity fluctuations is simi- 
lar in nucleon-nucleon interactions and central heavy ion 




10^ 
\|^ (GeV) 



FIG. 24: (Color online) UrQMD results of scaled variance uj 
of negatively charged hadrons in full phase space in inelas- 
tic p+p, p+n interactions and central Pb-fPb collisions as 
a function of collision energy compared to hadron resonance 
gas model predictions [lH for Pb-|-Pb collisions. The plot is 
taken from Ref. 



collisions. Thus with respect to the scaled variance of 
multiplicity distributions UrQMD behaves like a super- 
position model. The energy dependence of uj is different 
from the predictions of the hadron resonance gas model. 
UJ in UrQMD shows a strong increase with collision en- 
ergy in accordance to the experimental p+p data, while 
the hadron resonance gas model has a much weaker en- 
ergy dependence. 

In order to compare the UrQMD model to the exper- 
imental data, both the acceptance and the centrality se- 
lection of the NA49 experiment have to be taken into 
account. The predictions of the model, published in [l7|, 
are compared to the experimental data in Figs. [T^lfTil 

Two different centrality selections (see section IIII Bp 
are used in the model: first, collisions with zero impact 
parameter (open circles), second the 1% most central col- 
lisions selected in the same way as done in the experimen- 
tal data using a simulation of the acceptance of the veto 
calorimeter (full dots). 

The UrQMD model with collisions selected by their en- 
ergy in the veto calorimeter is mostly in agreement with 
data for all energies, acceptances and charges. UrQMD 
simulation of events with zero impact parameter (6 = 
0) gives similar results in the forward rapidity region, 
whereas uj is smaller in the midrapidity and the full ex- 
perimental regions, probably due to target participant 
fluctuations, which are still present for events selected 
by their forward going energy, but not for collisions with 
a zero impact parameter. 

The deviation of the multiplicity distribution from a 
Poisson distribution is similar in the model and in the 
data (see Fig. [^5]) . but the mean multiplicity is overpre- 
dicted in the UrQMD model for all rapidity intervals, 
charges and energies by about 20%. However, the scaled 
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variance of the multiplicity distribution is independent of 
mean multiplicity for superposition models. Since it was 
shown that UrQMD behaves like a superposition model 
for w, it is justified to compare lo for data and UrQMD 
even though the mean multiplicities are different. Within 
this framework one might speculate that the particle pro- 
duction sources in UrQMD are properly modeled but the 
number of sources is overestimated in central Pb+Pb col- 
lisions. 

In the experimental data an increase of fluctuations 
is observed when approaching midrapidity (Figs. [TBlllSp . 
The UrQMD model reproduces this behavior when a sim- 
ilar centrality selection is used as in the data. 

For the data an increase of u) is measured with decreas- 
ing transverse momentum at forward rapidity (Fig. I19p . 
In the UrQMD model a similar trend is observed, but lo is 
underpredicted at low transverse momenta. This might 
be related to effects like Coulomb and Bose-Einstein cor- 
relations, which are not implemented in the model. 

The HSD transport approach, following a similar strat- 
egy as the UrQMD model, yields similar results for uj. 
The energy dependence for central (6 — 0) Pb-f Pb colli- 
sions obtained by the HSD model are presented in [30| . 
These predictions were compared to preliminary NA49 
results on multiplicity fluctuations in j42| and were found 
to agree in the forward acceptance. Unfortunately HSD 
calculations for the larger acceptance used in this paper 
are not available yet. 



C. Onset of Deconfinement 



The fluctuations of thermalized energy are obtained by 
UrQMD and HSD simulations and are found to be 
5E/E < 0.03. 

Using this result one can estimate the additional 
multiplicity fluctuations of negatively charged hadrons 
caused by the fluctuations of thermalized energy to be 
^SE{h~) « 0.02 for the pure hadron gas or quark gluon 
plasma phase. In the mixed phase the expectation for 
^^SE{h~) amounts to ~ 0.03 at 80^ GeV. The predicted 
increase of lo by 0.01 due the mixed phase is smaller than 
the systematic error on the measurement of w. Therefore 
the data can neither support nor disprove the existence 
of a mixed phase at SPS energies. 

D. First Order Phase Transition 

It is suggested in [29ll that droplets of hadronic matter 
should be formed in matter when the system crosses the 
flrst order phase transition line during cool-down. These 
droplets are expected to produce multiplicity fluctuations 
10-100 times larger than the Poisson expectation in the 
full phase space. No predictions of the increase of uj for 
the limited experimental acceptance are available, but 
naively it can be expected to be of the order of 1-10 
(according to Eq. [3]). 

In our acceptance an excess of multiplicity fluctuations 
with respect to the UrQMD baseline, which does not in- 
clude an explicit phase transition, of larger than 0.1 can 
be excluded (see Fig. [25|) . 



In heavy ion collisions initial fluctuations in the 
stopped energy E are expected to cause fluctuations in 
the entropy S [l^. The energy dependences of various 
hadron production properties, like the kaon to pion ra- 
tio, the inverse slope parameter of kaons and the pion 
multiplicity [l^, |43| show anomalies at low SPS ener- 
gies which may be attributed to the onset of deconfine- 
ment [ll|. In [l2| it is predicted that this should lead to 
a non-monotonic behaviour of the ratio of fluctuations of 
entropy to stopped energy 



Rp — 



{SSf/S^ 
{5E)yE^ 



(11) 



At intermediate SPS energies, where a mixed phase of 
hadron gas and QGP is assumed, a "shark-fln" struc- 
ture with a maximum near 80A GeV is predicted. is 
approximately 0.6 both in the hadron and quark gluon 
plasma phase, in the mixed phase it can reach values up 
to 0.8. 

In [l^ these relative fluctuations are related to multi- 
plicity fluctuations under the assumption of a proportion- 
ality of entropy to produced particle multiplicity, namely: 



E^ 



■ in) ■ Re. 



(12) 



E. Critical Point 

It is expected that the hadron gas and quark gluon 
plasma regions in the phase diagram of strongly interact- 
ing matter are separated by a flrst order phase transition 
line at high baryo-chemical potentials and moderate tem- 
peratures. A crossover between both phases is predicted 
for high temperatures and low baryo-chemical potentials. 
Then the flrst order phase transition line will end in a 
critical point. 

If the freeze-out of matter happens near the critical 
point, large fluctuations, for instance in multiplicity and 
transverse momentum, are expected. In [Tsj it is esti- 
mated that the scaled variance of the distribution of total 
multiplicity of single charged hadrons should increase by 
about 1 near the critical point. However, this estimate 
has a large and difficult to estimate systematic error. The 
limited acceptance should reduce the critical point signal 
by a factor of about 2. Consequnetly, the expected in- 
crease of the scaled variance in the vicinity of the critical 
point is about 0.5. 

These critical fluctuations are expected to be located 
mainly at low transverse momenta [l3| . The scaled vari- 
ance as a function of the baryo-chemical potential is com- 
pared in Fig. [2n]to the UrQMD baseline. As the increase 
of fluctuations due to the freeze-out in the vicinity of 
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FIG. 25: (Color online) Top: Scaled variance lo of the multi- 
plicity distribution of negatively charged hadrons at forward 
rapidities produced in central Pb+Pb collisions as a function 
of the baryo-chemical potential /xs ,48] . In addition a sketch 
of the expected increase of u due to the critical point p^ . l49l ] 
is shown. The UrQMD results are given for a centrality se- 
lection similar to the experimental data. Bottom: Ratio of 
in data and UrQMD as a function of [i,b- 



the critical point is expected to be restricted to a range 
in the baryo-chemical potential which is comparable to 
the difference in baryo-chemical potentials of the different 
collision energies [50|, the signature of the critical point 
is expected to increase w at one collision energy only. A 
sketch of the expected increase of uj due to the critical 
point [13] is shown in Fig. [551 No significant increase 
of w which may be attributed to the critical point is ob- 
served in the data. The scaled variance for low transverse 
momentum particles (see Fig. [T5]) does not show a signifi- 
cant non-monotonic structure or excess over the UrQMD 
baseline either. 



158 A GeV was studied for positively, negatively and all 
charged hadrons. The total selected experimental ac- 
ceptance (0 < y(ji) < Ubeam) IS divided into a midra- 
pidity (0 < tj^tt) < 1) and a forward rapidity region 
(1 < y{TT) < ybeam)- At forward rapidity a suppres- 
sion of fluctuations compared to a Poisson distribution is 
observed for positively and negatively charged hadrons. 
At midrapidity and for all charged hadrons the fluctu- 
ations are higher. Furthermore the rapidity dependence 
at all energies and the transverse momentum dependence 
at 158^ GeV were studied. The scaled variance of the 
multiplicity distribution increases for decreasing rapidity 
and transverse momentum. 

The string-hadronic model UrQMD significantly over- 
predicts the mean multiplicities, but approximately re- 
produces the scaled variance of the multiplicity distribu- 
tions. 

Multiplicity fluctuations predicted by the grand- 
canonical and canonical formulations of the hadron reso- 
nance gas model [l^ overpredict fluctuations in the for- 
ward acceptance. The micro-canonical formulation pre- 
dicts smaller fluctuations and can qualitatively reproduce 
the increase of fluctuations for low rapidities and trans- 
verse momenta. However no quantitative calculation is 
available yet for the limited experimental acceptance. 

At RHIC and LHC energies the difference in w for the 
string-hadronic and the hadron-gas models in the full 
phase space is much larger than for SPS energies and 
experimental data should be able to distinguish between 
them rather easily. 

Narrower than Poissonian (a; < 1) multiplicity fluc- 
tuations are measured in the forward kinematic region 
(1 < 2/(7r) < ybeam)- They can be related to the reduced 
fluctuations predicted for relativistic gases with imposed 
conservation laws. This general feature of relativistic 
gases may be preserved also for some non-equilibrium 
systems as modeled by the string-hadronic approaches. 

The predicted maximum in fluctuations due to a first 
order phase transition from hadron resonance gas to 
QGP [l2| is smaller than the experimental errors of the 
present measurements and can therefore neither be con- 
firmed nor disproved. 

No sign of increased fluctuations as expected for a 
freeze-out near the critical point of strongly interac ting 
matter was observed. The future NA61 program 
will study both the energy and system size dependence 
of fluctuations with improved sensitivity in a systematic 
search for the critical point. 



APPENDIX A: DERIVATIONS 



1. Acceptance Dependence of u 

VI. SUMMARY 

Provided the particles are produced independently in 
The energy dependence of multiplicity fluctuations in momentum space and the form of the momentum distri- 
central Pb-f Pb collisions at 20^, 30^, 40^1, 80A and bution is independent of multiplicity, the scaled variance 
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in a limited acceptance is related to the scaled variance 
in full phase-space ("47r") by an analytic formula. 

Under these assumptions, having an experimental ac- 
ceptance registering the fraction p of the total number 
of tracks N is equivalent to roll a dice for each particle 
in the full phase space and to accept it with a probabil- 
ity of p. Therefore the probability to measure a number 
of particles n in a fixed acceptance follows a Binomial 
distribution: 



duction sources: 



B{n\N) 



Nl 



n!(7V-7i)! 



(Al) 



For a number of particles varying in the full phase space 
according to P47r (N) , the probability to measure a num- 
ber of particles n in the limited acceptance is: 



PA(n) = ^B(n|iV)P4.(iV). 



(A2) 



N 



From Eqs. I All A2l follow that the mean number of parti- 
cles in the acceptance is: 



< n >= p < N >, 



(A3) 



and the variance of the number n of particles in the ac- 
ceptance is given by: 

Var{n) =< Var{n\N) > +Var{< n\N >) 

=< Var{n\N) > +Var{pN) (A4) 
=< N > p{l -p)+ p^Var{N). 

Finally, the scaled variance in the limited acceptance uJacc 
is related to the scaled variance in the full phase space, 
ijJA-n, as: 



p {lU4t, - 1) + 1. 



(A5) 



The acceptance dependence given by Eq. IA5I is not 
valid when effects like resonance decays, quantum statis- 
tics and energy- momentum conservation introduce cor- 
relations in momentum space. 



2. Participant Fluctuations 

In a superposition model the multiplicity n is the sum 
of the number of particles produced by k particle pro- 



(A6) 



where the summation index i runs over the sources. As- 
suming statistically identical sources the mean multiplic- 
ity is: 



(n) = (k) (n-) 
and the variance reads: 



(A7) 



Var{n) = (fc) l/ar(n"°) -t- (n"")^ Var{k). (A8) 
Using these equations the scaled variance of n can be 
expressed as: 



1^ — , . , , — — : 1 ,. . , — — ; — U! -\- \n ) ■ UJk- 



(k) {n^") 



(fc) (n'' 



(A9) 

For the case of a constant number of sources the scaled 
variance is independent of the number of sources. 
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energy 




(GeV) 


< y(7r) < yh^am 


< j/(7r) < 1 


1 < J/(7r) < J/fceam 


20A 


0.88 ±0.02 ±0.02 


0.99 ±0.02 ±0.02 


0.86 ±0.02 ±0.02 


30A 


0.85 ±0.01 ±0.02 


0.96 ±0.02 ±0.02 


0.84 ±0.01 ±0.02 


4:(jA 


0.89 ±0.01 ±0.02 


1.01 ±0.01 ±0.02 


0.87 ±0.01 ±0.02 


oU/i 


0.93 ±0.03 ±0.02 


1.04 ±0.03 ±0.03 


0.89 ±0.03 ±0.02 


iOO-ri 


0.89 ±0.02 ±0.02 


1.00 ±0.02 ±0.02 


0.84 ±0.02 ±0.02 


energy 




(GeV) 


< ?/(7r) < ybeam 


< y(7r) < 1 


1 < J/(7r) < yteam 


20A 


0.94 ±0.02 ±0.02 


1.01 ±0.02 ±0.02 


0.93 ±0.02 ±0.02 


30A 


0.91 ±0.01 ±0.02 


1.01 ±0.02 ±0.02 


0.91 ±0.01 ±0.02 


4QA 


0.92 ±0.01 ±0.02 


1.02 ±0.01 ±0.02 


0.91 ±0.01 ±0.02 


an A 

oU/l 


0.88 ±0.03 ±0.02 


1.05 ±0.03 ±0.02 


0.86 ±0.03 ±0.02 


A 


0.90 ±0.02 ±0.02 


1.05 ±0.02 ±0.02 


0.83 ±0.02 ±0.01 


energy 




(GeV) 


< 1/(7r) < yhcam 


< j/(7r) < 1 


1 < J/(7r) < yteam 


20A 


1.01 ±0.02 ±0.04 


1.10 ±0.02 ±0.04 


0.94 ±0.02 ±0.04 




1.01 ±0.02 ±0.04 


1.07 ±0.02 ±0.04 


0.94 ±0.01 ±0.04 


40A 


1.10 ±0.01 ±0.04 


1.15 ±0.01 ±0.04 


1.01 ±0.01 ±0.04 


mA 


1.21 ±0.04 ±0.05 


1.22 ±0.04 ±0.05 


1.07 ±0.03 ±0.04 


158 A 


1.24 ±0.03 ±0.05 


1.20 ±0.02 ±0.05 


1.09 ±0.02 ±0.04 



TABLE V: Scaled variance of the multiplicity distribution of 
positively (top), negatively and all (bottom) charged hadrons 
as a function of energy. The first error is the statistical and 
the second error the systematical one. 




200 220 240 260 280 200 220 240 260 280 

n(h*) n(h*) 



FIG. 26: (Color online) Left: multiplicity distributions of pos- 
itively charged hadrons in full experimental acceptance in the 
1% most central Pb±Pb collisions from 2^A (top) to 158^1 
GeV (bottom). The dashed lines indicate Poisson distribu- 
tions with the same mean multiplicity as in data. Right: the 
ratio of the measured multiplicity distribution to the corre- 
sponding Poisson one. 
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FIG. 27: (Color online) Left: multiplicity distributions of neg- 
atively charged hadrons in full experimental acceptance in in 
the 1% most central Pb-|-Pb collisions from 20A (top) to 158A 
GeV (bottom). The dashed lines indicate Poisson distribu- 
tions with the same mean multiplicity as in data. Right: the 
ratio of the measured multiplicity distribution to the corre- 
sponding Poisson one. 



FIG. 28: (Color online) Left: multiplicity distributions of all 
charged hadrons in full experimental acceptance in the 1% 
most central Pb-|-Pb collisions from 20A (top) to 158A GeV 
(bottom). The dashed lines indicate Poisson distributions 
with the same mean multiplicity as in data. Right: the ratio 
of the measured multiplicity distribution to the corresponding 
Poisson one. 
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FIG. 29: (Color online) Left: multiplicity distributions of pos- 
itively charged hadrons in midrapidity acceptance in the 1% 
most central Pb+Pb collisions from 20A (top) to 158A GeV 
(bottom). The dashed lines indicate Poisson distributions 
with the same mean multiplicity as in data. Right: the ratio 
of the measured multiplicity distribution to the corresponding 
Poisson one. 



FIG. 30: (Color online) Left: multiplicity distributions of neg- 
atively charged hadrons in midrapidity acceptance in the 1% 
most central Pb-|-Pb collisions from 20A (top) to 158A GeV 
(bottom). The dashed lines indicate Poisson distributions 
with the same mean multiplicity as in data. Right: the ratio 
of the measured multiplicity distribution to the corresponding 
Poisson one. 
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FIG. 31: (Color online) Left: multiplicity distributions of all 
charged hadrons in midrapidity acceptance in the 1% most 
central Pb+Pb collisions from 20^1 (top) to 158A GeV (bot- 
tom). The dashed lines indicate Poisson distributions with 
the same mean multiplicity as in data. Right: the ratio of 
the measured multiplicity distribution to the corresponding 
Poisson one. 



FIG. 32: (Color online) Left: multiplicity distributions of pos- 
itively charged hadrons in forward acceptance in the 1% most 
central Pb-|-Pb collisions from 20A (top) to 158A GeV (bot- 
tom). The dashed lines indicate Poisson distributions with 
the same mean multiplicity as in data. Right: the ratio of 
the measured multiplicity distribution to the corresponding 
Poisson one. 
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FIG. 33: (Color online) Left: multiplicity distributions of neg- 
atively charged hadrons in forward acceptance in the 1% most 
central Pb+Pb collisions from 20^1 (top) to 158^1 GeV (bot- 
tom). The dashed lines indicate Poisson distributions with 
the same mean multiplicity as in data. Right: the ratio of 
the measured multiplicity distribution to the corresponding 
Poisson one. 



FIG. 34: (Color online) Left: multiplicity distributions of all 
charged hadrons in forward acceptance in the 1% most cen- 
tral Pb-I-Pb collisions from 20A (top) to 158^ GeV (bottom). 
The dashed lines indicate Poisson distributions with the same 
mean multiplicity as in data. Right: the ratio of the measured 
multiplicity distribution to the corresponding Poisson one. 



